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Abstract. The theorem studied is known before. Here is given a new proof. The 
proof has been part of a course material on Sobolev space theory with a special kind 
of outlook. The proof here is in accordance to this goal. At the same time several 
ideas of interest are shown that can be of general use. 



Introduction 

The definition of Sobolev space used here is for reasons to obvious later on that 
of the closure in the approprate Sobolev space norm of the set of functions with 
m times continuous derivatives on some open subset of iV-dimensional Euclidean 
space. 

The Weierstrass approximation theorem is stated below and concerns approx- 
imation of contionuous functions by polynomials. The higher order Weierstrass 
approximation theorem concerns approximation of m times continuously differe- 
tiable functions by polynomials in such a way that all derivatives up to order m 
are approximated simultanuously. The theorem is stated below. The theorem is 
previuosly known. Here is given a very interesting proof that makes use of several 
ideas. These and their interplay displayed is the real goal! That a straightforward 
proof has less complexity is beside the point here. 

The theorem holds for a cube as domain and with the open bounded domain 
denoted O in R^, then, if the domain allows the Sobolev functions considered to 
be extended to the full space, then the theorem follows for this domain as well. The 
argument is simply that take cube enclosing the domain and the do the extension of 
the function to full space and then restrict to the cube. Then use the approximation 
by polynom for the cube as in the theorem. Conclude by taking the restriction of 
the polynomial to the domain. 

A very good extension theorem is that of P.W. Jones [1]. The domains al- 
lowing extension by his method are called uniform domains or also by the name 
(e, 5)-domains. Domains with self-similar fractal boundary are simple examples of 
uniform domains. There is a point here that the Jones extension theorem also 
allows for p = oo as well. 

A consequence is that polynomials are dense in integer order Sobolev spaces with 
1 < p < oo when the domain is bounded and is a so called uniform domain. 
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Theorem and Proof 

Definition: Let O, be open in R^. For u a "nice" real function on O the Sobolev 
space norm is 

IMIiw(fi)= Yl II^NU^n) ~ (\\ Dau W P L*(n))* ■ 

\ce\<m \a\<m 

Definition: Let O, be open in R^. Define C m (0) as the m times continuously 
differentiable real-valued functions on O. 

Notation: We use V k u, the /c-gradient, in integral expressions as follows 

||V fc w|| L p ( n) = Yl IP^IU^n)- 

|a|=fc 

Notation: Let 7? be all polynomials in R^ and V k = {P e V : V k+1 P = 0}. 
Then some easy results without proofs. 

Theorem: A Poincare inequality, dim general, order one and general p. 

Let 1 < p < oo and Q = [0, 1] N . Let u e C°(Q) and D lU e C°(Q). Furthermore 
w|xi=o = 0. Then 

IMIlp(q) < \\D\u\\lv{q)- 
So, here Q is a unit cube and D\ is the derivative in the xi-direction. 

Corollary: The Standard Poincare Inequality. 

Let u G C™(Vt) for O bounded, open subset of R^ and \a\ < m. Let A = 

Y \\D a u\\ LPm < A\\V m u\\ LP{n) . 

\ct\<m 

Corollary: A more detailed Poincare inequality. 

Let Q = [0, 1]^ and let a and f3 be multi-indices that are partially ordered by 
"<" in the natural way. Let \a\ = t, \(3\ = m and a < (3. Assume that multi-indices 
{lk)T exist such that 7< = a and 7 m = (3 and also ^ t < It+i < • • • < lm-i < 7m- 

Let u E C m (Q) and 

(D^ k u\ x -o) = 

V fk-l~ U / 

or 

(D lk u\ x -i) = 0. 

Then for A = A Njrnjt , P ,Q, 

1 1 7-ia„. 1 1 /ill r->/3„.l I 
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Theorem: (Weierstrass Approximation Theorem.) 

Let Q C H N be connected, open and bounded. Then V is dense in L°°-norm in 



This will be a tool for proving a higher order version. 
The proof here is new. 

Theorem: A Higher Order Weierstrass Approximation theorem. 

Let m > 0, integer, and let Q = [0, 1]^. Then V is dense in VF m '°°-norm in 



The difficulty with the use of Weierstrass approximation theorem here is that 
there are many partial derivatives. They shall all be approximated by the partial 
derivatives of only one polynomial. - A trick (=method) is called for. 

First an observation. 

Notation: Let W m ' p (Cl) : where O is an open subset of R^, be the closure of 
C m (0) in the W m > p (Q)-norm. 

Theorem: Let Q = (0, 1) N then W m > p (Q) = W m > p (Q). 

Proof: Obviously W rn ' p (Q) C W m ' p (Q). Hence it is enough to prove the other 
inclusion. 

Let u G W m,p (Q). Then scale=dilate the variables with a factor 1 — e. Then 
u is transformed into u and u G W m ' p (Q i/i_ e ). Take restriction u\q. Clearly 
u\q G W rn ' p (Q). Let e take values {n -1 }, then {u n } a sequence which tends to u 
in the norm of W m,p (Q). 

End of Proof 

Proof: (Wannebo) The higher order Weierstrass approximation theorem. 

This proof can regarded as good training ground for the ideas discussed so far. 
It is not really meant for memorizing. 

The key to get started with the proof is to change the setup and study another 



We choose the space u G C Nm (Q) instead of C m {Q). 

Fix the special partial derivative D 13 = T[f =1 D™. This way the multiindex (3 is 
defined. 

Define the set S as the set of 2 Nm elements, each defined by how it acts on the 
monomial x 13 . All transformations, which to each one-degree factor in ', say the 
factor Xj, substitute with either Xj or 1 — Xj in this case. The set S is all possible 
combinations of such transformations. We write S = {a} 



C°(0). 



C m (Q). 



space. 



Identity: 
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- Check! 

Let u G C rnN (Q). Then by the Weierstrass approximation theorem, given e > 0, 
there is a polynomial Q a , such that 

\\DP[<j{xP)u] -Qa\\ L °°(Q) < e- 
Next step is to find a polynomial P„ such that 

\\DP[a{xP){u-P (J )]\\L~ { Q)<e. 
In order to prove this observe that the equation 

DP[<r(xP)P a ] - Q a = 

is solvable for P a for every Q a . The equation is linear and it is enough to solve for 
Per when Q a = x 1 any 7. Observe that 

D^[cr{x l3 )x 1 ] = const ■ x 1 + lower order terms. 

Hence it is possible by iteration to find P a . Solve with higher orders monomials 
first subtract, iterate. Since the equation is linear the general Q a has a P a solution. 

We are interested in the case \a\ < m and then it follows automatically that 
a < f3. This is the reason for studying C mN (Q) instead of C m (Q). 

Collect the results so far and use the Identity. Then 

\\D a [u-J2<r(^W\\L~(Q) = \\D a [J2<T(x P )(u-P*)]\\L~(Q) = 

a a 

= \\Y,D a [o-{x?){u - P ff )]|Uoo (0) < J2W D >( xf3 )( u - P °)]\\l~ {Q) < 
Use the more detailed Poincare inequality (the Corollary) 

< ^\\DP[a{ X P){u- P a )]\\ Loo{Q) < const. • e. 

This proves that V is dense in C mN (Q) in the VF m ' 00 ((5)-norm. 

In order to finish the proof, it only remains to prove that C mN (Q) is dense in 
C m (<5) in the IU m '°°(g)-norm. 

Hence let u G C m (Q), general. The space is complete in the VF m '°°((5)-norm. 
It is only needed now to find a Cauchy sequence 

{v n }?eC mN (Q) 

which converges in VF m '°°((5)-norm to u. As simplification, assume that Q has 
centre at the origin and side equal to 1. 
Make the coordinate transformation 

x — > (1 — l/(n + l))x, with n positive integer. We have the mapping Q — > Q n , 
where Q n has side 1 + 1/n. There is some convolution kernel with radial symmetry, 
\l/(r), which has support in a ball with radius 1 and with # G C m7V " m . Let 



A HIGHER ORDER WEIERSTRASS APPROXIMATION THEOREM - A NEW PROOF 5 



- a standard transformation. Let A n — > oo in fast enough. Now convolutions has 
such properties with respect to differentiation that 

Vn = ^n*U n \ Q eC mN (Q). 

But C m (Q) gives uniform continuity for the derivatives up to order m so it follows 
(Check!) that {v n } is a Cauchy sequence with limit u. 

End of Proof. 
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